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Instanton analysis is applied to model A of critical dynamics. It is shown that the 
^O ' static instanton of the massless 0^ model determines the large-order asymptotes of 
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^ , 1 Introduction 



The knowledge of large-order asymptotic beliaviour of perturbation series of 
static field-theoretic models is important for resummation of series for criti- 
cal exponents and scaling functions [1]. This behaviour has been thoroughly 
explored with the aid of instanton analysis and applied to the resummation 
problem in the prototypical static cj)^ model [2,3,4], which has been widely 
used as a model of critical behaviour in continuous phase transitions of ferro- 
magnetic type. 
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However, little is known about large-order asymptotics in dynamic field the- 
ories constructed from Langevin equations with the aid of the Martin-Siggia- 
Rose (MSR) formalism [5]. Recently, it was stated [6,7] that there is no instan- 
ton within the MSR approach in the Kraichnan model - which has attracted 
considerable attention as a model describing intermittency in turbulent diffu- 
sion [8] - and that the method of steepest descent has to be used in Lagrangean 
variables. 

In this paper we propose a method to assess large-order asymptotic behaviour 
of dynamic models near equilibrium, i.e. with Gibbsian static limit (we also 
will restrict ourselves to models generated by Langevin equations without 
mode coupling). We will use the steepest descent method (instanton approach) 
to find the large-order behaviour in a representative model. We will discuss 
large-order asymptotes of correlation functions, response functions and critical 
exponents. For equal-time correlation functions this asymptotic behaviour is 
shown to coincide with that of the static instanton approach. 

The particular model we deal with in this paper is one of the standard dynamic 
0^-based models: model A in the classification of Ref. [9]. In this model critical 
exponents are the same as in the static 0^ model apart from the dynamic 
exponent z. The large-order asymptotics of the dynamic exponent have not 
been analyzed so far. It should be noted here that the use of Lagrangean 
variables becomes prohibitively difficult in this case due to the essential non- 
linearity of the problem (the Kraichnan model is linear in the advected scalar 
field). 

The present article is organized as follows: in Sec. 2 construction of the MSR 
field theory corresponding to a non-linear Langevin equation is reviewed with 
special emphasis on the treatment of the functional determinant for the steep- 
est descent method. Existence of the dynamic instanton and relation of the 
dynamic instanton solution to the static one is analysed in Sec. 3. The fiuc- 
tuation determinant is calculated in Sec. 4, whereas Sec. 5 is devoted to a 
brief analysis of correlation and response functions. Results of this paper are 
summarized in Sec. 6. 



2 Dynamic field theory 



Consider the Langevin equation 



— ^ + r — = ^. 

dt 6(p 



where ^ is a Gaussian random field with zero mean and the correlation function 

(e(t, x)e(t', x')) = D{x - x') = 2T6{t - t')5(x - x') , 

where the shorthand notation x = (t, x) has been used. In equation (1) the 
action is the static action of arbitrary model with the known instanton. The 
most interesting example is the massless (^^ model: 

s = \d^d^ + 1^^ . (2) 

Here and henceforth, all necessary integrals and sums are implied. We find it 
convenient to start the analysis on a finite time interval [to,T], but in the full 
Euclidian space R^. In the MSR approach, functional integrals for correlation 
and response functions are calculated with the "measure" 

m[i^, v?'] = V^V^' det M e^^+^^'^+V'^' , (3) 

where the dynamic action is of the De Dominicis-Janssen form [10] 

The operator determinant det M in (3) has different interpretations. Most 
often it is written as a result of a formal change of variables in a 5 functional 
imposing condition (1), which yields [1,10] 

_ d5{x-x') ^ 5^S 



dt 5(p{x)6(p{x') 

for the operator M, after which the determinant is calculated in a formal loop 
expansion as 



det M = det I — + Tk]6{x- x') exp - f f dxdx' 



6'Sr 



5(p{x)5(p{x') 



(5) 



where K is the differential operator of the free-field part of the static ac- 
tion [e.g. K = — V^ for action (2)] and Sj its interaction part. An annoying 
feature of this straightforward method is that the normalization determinant 
det [{dt + rK)6{x — x')] = detMo is that of a non-self-adjoint operator. The 
explicit content of this factor is unimportant in the ordinary perturbation the- 
ory, because in the calculation of all expectation values it simply cancels out in 
each term separately. However, in instanton calculus the saddle-point solution 
which serves as the expansion point for a steepest-descent calculation of the 
functional integral is different from that of the ordinary perturbation expan- 
sion. Therefore, the possible cancellation of determinants requires additional 
analysis, which is much more convenient to carry out in terms of determinants 
of self-adjoint operators. 



It is quite possible and completely consistent to arrive at a functional measure, 
in which the normalizing determinant is that of a self-adjoint operator. To this 
end, recall the basics of the MSR construction, in which classical quantities are 
replaced by field operators '0, il)' with the usual bosonic commutation relations 
{normalized as [■?/'(t,x),-?/''(t,x')] = 5(x — x')} and the dynamics are given by 
the operator equations [5] 



dil> 
'dt 



where the MSR " Hamiltonian" is 



dip' 
'dt' 



ij',n 



n = --^'D^' + Tij' 



6ip 



(6) 



ip^-i/) 



We have chosen the MSR Hamiltonian in the form (6), which is slightly differ- 
ent from the originally proposed (although this possibility was discussed in the 
original paper [5]), because we want to write the free- field part of the dynamic 
action in the form of a formally convergent Gaussian integral at the outset. 
The chosen MSR Hamiltonian corresponds to the stage, where all integrations 
over the random source C, of Langevin equation (1) have been carried out, and 
- integrated over time - it yields the dynamic action (4). 

With the use of the traditional operator formalism the generating functional 
of Green functions of these fields may be expressed in the functional form [11] 



G{A^,A^>) = exp 



-—A— 
25¥ 5¥ 



exp 



Si + A^iP + A^,iIj' 



•^=0 



(7) 



where Sj is the time integral of the interaction part of the MSR Hamiltonian 
(6) in the interaction picture, whereas the kernel of the functional differential 
operator - written in terms of the two-component field \E' = {ip,ip') - is the 
propagator 



A{x,x') 



Aii(a;, x') Ai2(a;, x') 
A2i(x,a;') 

{dt + TKy^ D {dt + TK)-^] {x, x') {dt + TK)'^ {x, x') 
{dt + TK)-^{x,x') 

Mq^DMq^] {x,x') Mq\x,x') 
Mo^{x,x') 



In this approach the functional integral arises from expressing the functional 



differential operator in (7) as the Gaussian functional integral: 



exp 



-—A— 
2M M 



det M^Mq 



X [V(pV(p'exp 



.l*» + ^|; + ^'A 



. (9) 



where $ = {^p, (/?') and the operator /C satisfies /CA = 1 with the appropriate 
boundary conditions. This is the only place where a functional determinant 
arises and now it is that of a self-adjoint operator MqMq completely deter- 
mined by the free-field part of the MSR Hamiltonian. 

It might be asked what happened to the field-dependent extra factor in (5)? 
The answer depends on the definition of the time-ordered product (T prod- 
uct) at coinciding time arguments. If the time-ordered product at equal times 
is defined as the symmetric product of the field operators, then the inter- 
action functional Si in (7) is the Sym form [11] of the interaction part of 
the MSR Hamiltonian. The Sym form of the MSR Hamiltonian (6) contains 
terms quadratic in fields which exactly correspond to the field-dependent fac- 
tor in (5). It is quite possible, and here in fact more convenient, to define the 
coinciding-time time-ordered product as the normal product [11] of the field 
operators. In this case the interaction functional Si is taken in the normal 
form, which in the perturbation theory implies the condition that there are 
no graphs with closed loops of the retarded propagator {dt + TK)~^{x, x') at- 
tached to the interaction vertex. This amounts to vanishing retarded free-field 
propagator at coinciding times, and no additional factors to the determinant 
detMj^Mo are generated. Henceforth, we will use the latter approach which 
allows to choose the determinant factor of (3) in the field-independent form 

det M = J det M^Mq. 



In the instanton approach the space of integration in the functional integral 
is chosen using the properties of the saddle point and may well be different 
from that of the ordinary perturbation theory. Since the starting point is the 
functional integral constructed for perturbation theory, let us recall what the 
appropriate space of integration is in that case. In the proper choice of the 
formal space of integration there are two main points to be taken into account 
[11]: the operator /C must be nondegenerate and the kernel of the uniquely 
defined operator /C~^ must be the propagator A used for the construction of 
the perturbation theory. For the model generated by the static action (2) the 
usual diffusion kernel in unbounded space 



dt - rv^) (x, x' 



e{t- t') exp 



.'\2 



'4r(t-t') 



[4vrr(t-t')] 



d/2 



is the choice to generate the propagator A in (8). The space of integration 
E{A) conforming to these requirements may be constructed as 






An A 
A21 




An^ + A12V 

A21^ 



(10) 



where 77 and rj' are fields vanishing at the boundaries (i.e. at |x| -^ 00, t = tg 
and t = T) to ensure all partial integrations in the differential operators of /C 
without surface terms. It is worth noting that ^p'{T) = 0, but both fields have 
finite (and real- valued) initial values at t = to. 

In terms of fields vanishing at the boundaries the right-hand side of equation 
(9) assumes a bit different form giving rise to the representation 



exp 



- — A — 



det A12 A21 / Vrj'Vrj exp 



'-r]Aiir] - r]Ai2r]' 



X exp 



(Anr, + A,W) ^ + ^2iVj^ 



:iii 



in which the Gaussian functional integral is taken over fields vanishing at 
boundaries. Nevertheless, formal convergence requires it to be understood as 
an iterated integral in the order indicated in (11), the field r) to be real and 
the field 77' purely imaginary. 



3 Instanton analysis 



Let us calculate the parametric integral expressing the A^th order contribution 
to perturbation expansion in g of the A;-point correlation function 



1 

27[i J g 



I?^P(^V(xi, ti) . . . v^(xk, h) e-^ - ^Ig^ 



ffVifVif' ( 



-On 



(12) 



by the method of steepest descent in the variables g and (/?. We have written 
the normalizing determinant in the form of the functional integral 

(detMo^Mo)^^^^ = IIVipVip'e~'^o 
with the free-field part of the dynamic action 



in order to make explicit the cancellation of Jacobians accompanying several 
forthcoming changes of variables. We will show that the instanton solution for 
expression (12) may be constructed in close relation to the instanton solution 
of the corresponding equilibrium model. Therefore, it is worth reminding that 
the equilibrium expression for the A;-point correlation function without vacuum 
loops [there are no vacuum loops in the dynamic correlation function (12)] of 
the field y? is recovered in the limit to ^ —oo for ti = t2 = ■ ■ ■ = tk- 

Applying the instanton approach we make - in analogy with the usual instan- 
ton analysis of the ip^ model - the change of variables if -^ viV(/9, cp' -^ yNip\ 
g -^ g/N; in case of any other model it is convenient to scale the field variables 
extracting the dependence of N in the expression for S. The corresponding 
Jacobians from the numerator and the denominator of (12) cancel. The sta- 
tionarity equations for the exponential of the numerator are 



^6^S 



,v=o, 



5S^ d^' 

Sip dt ' 6ip' 

SS ^ dip 5S 

as __i 

dg g ' 



(13) 
(14) 
(15) 



To construct the instanton solution, consider the following auxiliary equation 
[note the change of sign of the time derivative in comparison with the the 
Langevin equation (1) and the corresponding stationarity equation (14)] 



dt 6(p 







(16) 



and assume for the time being that there is a solution ipd conforming the 
Cauchy condition for the finaltime instant ipd{T, x) = v^/(x) with a sufficiently 
smooth function v^/(x). Substitution of this solution to equation (14) yields 






-Dip' + 2 



d(pd 
dt 







v^vd 



which leads to the nontrivial 



V'd = ^D 



-I dip a 



dt 



(17) 



On the other hand, equation (13) on this solution is satisfied automatically: 
indeed, taking the time derivative of equation (16) yields 



ddipd_^_5'S 



dt dt 



5ip5(p 



d(pd 

dt 

""Pd 




v-^Vd/ 



difd 
dt 



which, apart from a multiphcative factor, is equation (13) on the solution (17). 

Finally, equation (15) imposes an asymptotic condition on ipd- it may be cast 
in the form (we remind that on the left-hand side of this equation integrations 
over time and space in the action are implied and intact) 



as _ , d_5S_ 



^ d — — 
dt dg 

Lp-^ip^ to 



¥'— x/'d 



t^T dg *-**() dg g 



since Tip'^ = 2TD~^dt^d = dt^Pd- Taking into account the boundary conditions 
imposed on the solution of equation (16) we immediately see that on the 
solution (16) and (17) the third saddle-point equation (15) reduces to the 
following equation for the limiting values of the dynamic instanton 

dg dg g' 

of which ipf is the condition for the Cauchy problem and v^o(x) = ipd{to,x). 
Since equation (16) gives rise to iterative solution with the advanced diffusion 
kernel, we see that sufficiently small "initial" field (fo vanishes in the limit 
to -^ — C)0 and we thus recover the usual equation for the static instanton 

dS(ipf) 1 , , 

dg g 

which allows to identify the so far unspecified Cauchy value of the dynamic 
instanton with the static one: ipd{T,x.) = v^/(x) = (pst{x.). 

Substitution of the solution ipd and ip'^ determined by equations (16), (17) and 
(18) in dynamic action (4) leads to the remarkable result that the dynamic 
action on the dynamic instanton solution asymptotically coincides with the 
static action on the static instanton: 

S{Vd, v'd) = S{Vst) - S{^o) > S{^st) ■ 

to — ^— oo 

Thus, we arrive at the conclusion that the exponential factor in our steepest 
descent analysis of correlation function (12) as well as pre-exponential factor 
asymptotically are the same as in the corresponding equilibrium static theory, 
the difference being generated by the (fixed, no averaging assumed) initial 
condition for the field (fo only. 

Let us explain the choice of the auxiliary equation (16). To consider the asymp- 
totic behaviour of the instanton (y?^ at t — ;► ±oo we have to study dynamic 
properties of this equation by investigation of its fixed points. Apart from 



the assumed solution Lp^ discussed so far, equation (16) has two rather ob- 
vious time-independent solutions: the trivial solution (/? = and the static 
instanton (^gt which obeys dt'^at = and {5S/5ip)ip^^^^ = 0. The trivial solu- 
tion does not conform to the non-vanishing initial condition, but gives rise to 
T5'^S/5ip^ = TK. This is a positive definite operator which renders the trivial 
solution stable in the stationary limit to -^ — oo. The second variation of the 

dynamic action for the stationary instanton is F I i^ — -^V^^i ) • This operator 

is well known from the instanton analysis of the usual ip'^ model [1]. It has at 
least one negative eigenvalue which determines the direction along which the 
solution approaches the stationary instanton at t — > oo. Such a situation with 
the time-independent solutions is a heavy argument in favor of the existence 
of the dynamic instanton which behaves as 

lim ipd = ^0 ' , lim (pd = ifst , 

t — Hq t{) — *— oo t — >oo 

even when the explicit solution is not known. 

As an example of this statement, consider a simple dynamic model, viz. the 
zero-dimensional (p'^ model. The dynamic action corresponding to Eq. (1) for 
the "field" (p{t) (without x dependence) is S* = Sy?^ + 4^^^- The equilibrium 
limit of such a model is the famous one-dimensional integral 



/ dxexpl— — a; — jr x 



which is often considered as an example of applicability of the instanton ap- 
proach to the large-order asymptotic investigations [1,2]. Then equation (16) 
assumes the form 



dtip-T[Tip + ^ip^]=0 



with the solution 



^' ^° V (6r + gy^l) exp [-2rr(t - to)] - g^l ' ^^^^ 



This dynamic instanton tends to (pst = J^T~/{—g) at t — ;► oo and vanishes at 

t -^ — oo. 

As to the usual ip'^ model in 4 — e dimension, we have not been able to find a 
closed-form representation for the dynamic instanton, but a convergent iter- 
ative solution of equation (16) may be constructed as follows. As the zeroth- 
order approximation, take 



V^?^(^.x) = jd^ 



exp 



(x - xQ^ 
'AT{T-t) 



[4vrr(T - 1)] 



d/2 



Vst{y^ 



which obviously is a solution of the linearized equation with the final condition 

(f^ {T,x.) = V9st(x). Rewriting equation (16) as the integral equation 



V9rf(t,x) = — ^|rft'|rfx' 



exp 



.'A 2 



[X — X. 

'4r{t' -t) 



^v^2(t',x)+v^r(t,x 



(0), 



[47rr(t' - t)f^ 
we immediately see that the next-to-leading term of the iterative sequence is 



T exp 

t 



J\2 



(^X — X ^ 

'AT{t'-'t) 



[47rr(t' - t)] 



d/2 



M 



v7{t\^) 



with the vanishing final value: tp]i {T,x) = 0. Proceeding in the same fashion 
we arrive at the usual tree-graph solution of the nonlinear equation (16) with 
the given final Cauchy condition instead of the initial condition. 

Similar solutions were used in [13,14] for the massive ip'^ model in dimension 
4 — e in the static limit. There the tree-graph expansion was sufficient to 
determine the asymptotic behaviour of the instanton at both large and small 
values of the coordinate argument |x|, which were of practical interest. 

It should also be noted that - unlike the usual field-theoretic perturbation 
expansion - the number of the graphs at large order of the tree expansion has 
no factorial growth. Therefore, this expansion is convergent in a disk of finite 
radius and can be analytically continued. For illustration of this statement, 
consider the expansion in the parameter g of the exact solution of the toy 
model (19), which obviously has a finite radius of convergence. 

Usually we are interested in the logarithmic theory. In this case the model is 
translation invariant with respect to x and scale invariant with respect to a 
synchronous dilatation of x and t (the scaling dimension of t being twice the 
scaling dimension of x). Therefore, the instanton has an arbitrariness of the 
form 



V^d(x,t) 



,D/2- 



-fd\ 



x 



Xo 



t 



(20) 



where D is the space dimension, xq and y are arbitrary parameters. The coeffi- 
cient r in (20) is introduced to make the second argument of fd dimensionless. 
The translation invariance in space and dilatation invariance forces to lift the 
corresponding degeneracy in the instanton solution. Due to the non-vanishing 
value v^^(T) = 2D~^dt(fd\t=T of the auxiliary field of the dynamic instanton 
solution, a shift of the integration space from the original perturbative inte- 
gration space (10) with <f'{T) = to a space with ^p'{T) 7^ is implied in the 
functional integrals of the /c-point correlation functions (12). 



10 



In this integration space the degeneracies of the instanton may be removed 
with the aid of the same unit decomposition 



1 = /rf^xo / d\ny'6 



xS 



D 



9_ 

24 



^ /c?xv9^(x)ln 



yrfxv9^(x)(x-Xo 



/ 


X-Xo\ 




[ y J \ 




) 


24 jd^vtoi^) 



D+l 



f211 



as in the static instanton theory, which then imposes conditions on the fi- 
nal value of the integration field: (/^(T, x) = (/9oo(x). {Relation (21) is written 
for a scalar field ip, in a vector case one more contribution due to the rota- 
tional invariance is necessary [1,13]}. The S functions introduced by the unit 
decomposition (21) fix the arbitrariness of the instanton completely. 

An appropriate change of variables to avoid the spatial translation invariance 
and dilatation invariance problem is well known [1,2,13]. Therefore, we will 
henceforth consider only specific dynamic properties of the model. All needed 
operations to deal with the spatial translation invariance and dilatation in- 
variance problems, together with spatial arguments of of fields and spatial 
integration, are implied and omitted in the following. 

Thus, we arrive at the correlation function 

1 ^ ^g //P(pPv^V(xi, ti) . . . v^(xk, 4) I e-^ - ^ Ig^ 



27fi J g 



ffv^V<p'e~^o 



where I stands for all contributions from the unit decomposition (21). At the 
leading order in N we replace all v^(xi, tj) by ipd in the pre-exponential factor. 
In the stationary limit to -^ — oo the initial value of the instanton ipd{to) 
vanishes and in the exponential the static instanton action is recovered at 
leading order in A^. 



4 Calculation of the fluctuation determinant 



During the calculation of the fiuctuation integral the usual stretching of vari- 
ables 6ip -^ 6ip/y/N, 5ip' -^ 5ip' /y/N in the fiuctuation integration in the nu- 
merator is accompanied by a similar transformation cp -^ Lp/yN ^ ip' -^ ip' /yN 
in the denominator: the Jacobians thus cancel. We shall not write the factor 
due to fiuctuations of g explicitly, because after the standard change of vari- 
ables in the original dynamic action ip -^ v/ \/—g-, V^' -^ v' / \/~g fiuctuation 
contributions due to fields 5ip and 5(p>' on one hand and due to the couphng 



11 



constant 5g on the other factorize, the latter giving rise to a trivial numerical 
factor. 

The fluctuation integral over 5ip and 5ip' may be written as 



(5S = J ffvS^VS^'exp 

X llv5ipV6ip'lexp 



^6^'D6^' - 6^' {^ + VK5^ 



1 1 8^S 

-5ip'D5ip' - -5ipTip'a-^ 



6ip 






5ip 



m) 



For economy of notation, we have used the same symbols for the integration 
variables 5^p' and 5^p in both the denominator and the numerator of (22). On 
the dynamic instanton we immediately see that 



5^S 



Sip' 



Wd 



6^S 



Vd 



5ip' 






fd 



With the aid of this relation and the change of variables 



^^ = ^' + 2D" 



dt 



(23) 



the quadratic form of the exponential of the numerator of fluctuation integral 
(22) may be expressed as 



\ €.'S = I ^'D^' - ^' (-i + rK + r^^ 



ip,ip' 



dt 



1 6'S 



6ip 
fd 
1 

0(Pst + 



<fst 






6^0, (24) 



Vo 



where 6ipst = Sip\t=T and 6ipo = 6ip\t=to- Note the change of sign of the 
time derivative in the second term of (24). We also remind that the term 



■5ip, 



St 



S'S 



5ip' 



6ipst in (24) describing the effect of asymptotic (equilibrium) fluc- 



fst 



tuations would be absent in the original perturbative integration space. 

We intend to calculate the fluctuation integral over the dynamic fluctuations 
regarding the term ip'T / 



6^' 
action, as perturbation. At 



6ip, where Si is the interaction part of the static 

fd 

eading order in A^, when only the fluctuation de- 
terminant is needed without any other fluctuation contributions, this trick 
leads to a sum of closed loops of advanced diffusion kernels which essentially 
reduces the purely dynamic fluctuation determinant to a constant regardless 



12 



of the properties of the full differential operator in (24). The account of dy- 
namic fluctuations becomes highly non-trivial, of course, when pre-exponential 
fluctuation terms are included. 

To arrive at a tractable perturbative expansion of the dynamic fluctuation 
determinant we write the normalization factor in terms of the same perturba- 
tion expansion as indicated in (24). Change of variables (23) casts the free-field 
action in the denominator in the form 



1 

2 



.lSl,^,So = ^^'D^'- 



,p,if' 



1 
2 



-■ip' f-^ + TK\ 6ip--6(f,tK6ifst + ^SifoK6ifo (25) 



with the same free- field part for the perturbative calculation as in (24). A shift 
of the integration space similar to that in the numerator of the fluctuation 
integral (22) is implied in (25) as well. 

The integration space corresponding to perturbative expansion with (25) as 
the free-field action includes fields constructed with the aid of the advanced 
diffusion kernel by the standard prescription, i.e. 



S^\ _( ^11 ^21 

V-' j I Ai2 




Anr] + A21V 



(26) 



where 77 and rj' are fields vanishing at the boundaries (i.e. at |x| -^ 00, t = T 
and t = to). However, the propagator structure leads to final values of the 
fields 6(f and ip', which become arbitrarily small, when T grows. Since in the 
integration space around the instanton the long-time asymptotics of the fluc- 
tuation field 6ip — * Sipst is finite and should eventually be independent of the 
choice of T, an additional contribution is required to account for the asymp- 
totic fiuctuations described by 6ipT- The spatial convolution of the advanced 
diffusion kernel with the difference of the desired asymptotic field 6ipst and 
the final value of S(f will do: 



6(pT{t,x.) = / dx.' 



exp 



X — X 



/\2 



4r(T - 1) 



[47rr(T - t)f^ 
since it is the solution of the homogeneous equation 



5^sM') - 6v{T,x'] 



(27) 



-^^(-^^ 



6ipT{t,x.) = 



with the final value (5(/9r(7', x) = (5(/9st(x) — 6ip{T,x.). 

Due to this, on the finite interval [to, T] the free- field action on the integration 



13 



field composed as a sum of (26) and (27), i.e. 

5ip = 5ip + 6(pT = All?] + A21V + 6(pT 
amounts to 
1 



(2J 



~ 2^'p,'p'^o 



l^'D^p'-ij'L^^+TKU^ 



- ^ 6MT)K6MT) + \ [5^0 + 5Mto)]K[5^o + ^Mh)] ■ (29) 
Here, liowever, lim SipxiT) = Sipst by construction (27). Tlierefore in (29) 

T — >oo 

tlie explicit dependence (apart from the upper limit of all time integrals) on T 
remains only in the initial value of the field 6(fT{to), which is finite. Moreover, 
from (27) it also follows that 6ipT{to) -^ 0, when T -^ 00. Thus, we may write 
the original quadratic form on the ray [to, 00) as 



2^'P,'P'^^ 



mn j i *W -*'(-! + FA-) *^ 



- - 5(pstK5(pst + - 6ipqK6ipq 

- an expression which in the functional integral of the denominator of (22) 
allows to factorize the asymptotic fluctuations described by 6ipst from the 
"purely dynamic" fluctuations described by 6ip and ip' from (26). 

A similar representation in the quadratic form of fluctuations (24) leaves ex- 
plicit dependence on 6lpt also in the "interaction" term: 



-I ^l.'S = ,iim^| I ^'Di^'-i'' (-| + ri^) ^-^' f r^^ 



-^^^(^)^ 



dipriT) + -difrit^ 



5ip 



6ip + 6ip' 



•T 



Vd 



ifist 



V 



Sipritc 



v>o 



Again lim 6(fT{T) = Sipst and lim v^t(^o) = by construction. Moreover, as 

T — ^00 T — ^00 

seen from (27), the latter limit applies for almost all t as well: lim frit, x)) = 

T — >oo 

0, to <t < T. From this it follows that 



^' T 






6(pT 



^'(t) r 



to 






1 — ^oo 

fd 



Therefore, we see that representation (27) and (28) leads to factorization of the 
time-independent asymptotic fluctuations Sipst and the dynamic fluctuations 
Sip and ip' in the numerator of (22) as well. 
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Therefore, the fluctuation integral may be expressed in the explicitly factorized 
form 



5S 



6'S 



X UvI^Vip' exp 



-i,'Di,'-i,' 



d5^ \ 1 




-ip'Dilj' - V' 



d 



dt 



TK\6^ 









(30) 



with the dynamic fields of the structure (26). On the right-hand-side of (30) 
the first line yields the familiar static fluctuation determinant including the 
unit decomposition (21). 

A non-trivial dependence on the initial time instant to remains in (30) through 
the non-zero initial value of the fluctuation field SipQ. This is an integral of 
the independent variables r] and r]' with vanishing boundary conditions from 
representation (26): 

oo 

6(Pq{x.) = dt Mx [Aii(to,t,x-x')?7(t,x') + A2i(to - t,x - x')?7'(t,x )] 
to 

and thus leads to fairly complicated propagator structure of the perturba- 
tive expansion of the dynamic contribution in (30). This also means that a 
functional dependence on the dynamic instanton remains in the fluctuation 
determinant for finite tg- In the asymptotic limit to -^ — oo - which restores 
the translation invariance with respect to time - the dependence on the in- 
stanton vanishes, however. This takes place, because due to the properties of 
the diffusion kernel, in the limit to ^ — oo the initial value of the fluctuation 
field vanishes: Sip^ -^ 0. Therefore, in this limit the dynamic part of (30) may 
be expressed in a form similar to (7): 






VS^Vip'exp 
ffvI^V^lj'lexp 



l^'D^p'-ij' 



■1f-H 



i^.W -./(-! + FA-) J^-V/(r0' 



•fid 



exp 



1 ^ AT ^ 

A^ — 

2(5^ (5^ 



exp 



-^' r 



6^' 



'Pd . 



■ (31) 



'I'=0 



Here A^ is obtained from the propagator kernel (8) by replacing the retarded 
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diffusion kernel by the retarded one and vice versa. Inspection reveals that 
the perturbation expansion of (31) contains only closed loops of the advanced 
propagator which all vanish (including the closed single loop, due to our def- 
inition of the T product at coinciding times). Therefore, we immediately see 
that perturbatively the expression in (31) is equal to unity and the dynamic 
fluctuation determinant (22) is equal to the static one in the limit to -^ — oo: 



V6(pstl exp -2^^/^3*7-72 



lim (5S — , 

to^ — 00 f f L 



1. 6^ 

6ip 




/D.,.exp(-i.,„A-.,„ 



(32) 



Thus, not only the pre-exponential and exponential factors in the large N 
asymptote of (12) coincide with the static expression, but the fluctuation 
integrals with the chosen regularization of determinants coincide as well. 



5 Dynamic part of the correlation function and response function 



We have shown that dynamic correlation functions coincide with the static 
ones at leading order in 1/A^ in the limit to ~^ —00, when the model is trans- 
lational invariant in time. Leading-order contributions to response functions 
vanish in this case, however. To see this, consider the A^th-order contribution 
to the response function 

1 , //l?V9l?V9V(xi,ti)v9'(x2,t2)Ie-^-^lg^ 
— I^IL = . (33) 



ffVifVip' e'^o 



On the instanton solution (16) and (17) the pre-exponential factor in (33) 
includes the factor 

d 

<^d(x2,t2) = 2D-^—<^d{^2,t2) 

which vanishes in the limit to -^ —00 due to the necessary asymptotic prop- 
erties of the instanton, which approaches a finite time-independent limit at 
long times after the initial time instant. Thus, the leading contribution to the 
response function is given by corrections of order 1/N. 

To capture dynamic effects it is sufficient to consider only corrections effected 
by the replacement ip ^ 6ip and ip' -^ S(p'. 1/N corrections of any other type 
lead only to the small changes to the stationary results obtained previously 
(without any dynamic dependence). 
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The fluctuation integral for the pair correlation function is 



X VScpVSif' 6(p{:x.i,ti)6(p{'K2,t2) I exp 



^s^'DSv' - 5^' (^ + ri^^v.) 



-1 



-V I — + r- 



dt 5ip' 



Sip 

•fid 

6ip 

•fid . 



With the use of the change of variables (23) and the subsequent decoupling 
of dynamic and static fluctuations according to (26), (27) and (28), this fluc- 
tuation integral gives rise to the sum 



VS(pstS(pst{^i)S(pst{^2) I exp "2 V 



SE 



6^ 



Sfst 

fst 



ipip 



+ 



VScpstew i--S(pstK6ipst 
IvSipstexp (--6ipstK6ipstj IvSipsd exp --S^p. 



6'S 



'■'W 



SiPst 

fst 



IJVdipVi)' exp 



l^'D^l,' - i,' L^ + VKbi^ + i b^^Kbip^ 



-1 



X V5(pVip' 6(p{:x.i,ti)5(p{:s.2,t2) exp 



^^'D^'-^'L^^+Tk\6^ 



-'(#: 



6(p+ - 6(p, 



6^S 



^d 



°V 



(5v9o 



Vo 



where the flrst term is the static fluctuation correlation function and in the 
second term the dynamic part - multiplied by the static fluctuation determi- 
nant (32) - may be expressed formally as 



A 



12 



1 + r 



5''Si 



dip' 



Ai2 

•fd 



D 



l + AoiT 



5''Si 



dip' 



H>d, 



^21 



in the limit t, 



-^ -oo- 



6 Conclusion 



In this paper flrst steps have been taken towards large-order (A^ -^ oo) asymp- 
totic analysis of non-linear stochastic fleld theories with the aid of the instan- 
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ton method. It is shown that this approach can be used in dynamic models as 
well as in the static ones. 

For the near-equilibrium model A we have shown that in the temporal translation- 
invariant case asymptotic properties at leading order in A^ are almost com- 
pletely determined by the corresponding static instanton solution and fluctu- 
ations thereabout which leads to factorial behaviour of the large-order con- 
tributions similar to that in the static instanton analysis: viz. the large-order 
asymptote for an arbitrary quantity F behaves as 

p[N] ^ N\Ca^N\ (34) 

where F^^^ is the A^th order of the expansion of F in the parameter e (e 
is the coupling constant g or the dimensional regularization parameter e), a 
and h are constants known from the analysis of the static theory [1,2], C is 
an amplitude factor - either a constant or a function of coordinate and time 
arguments depending on the quantity F . In the static (p^ model the functions 
C for universal quantities have been calculated in Ref. [14]. 

In our treatment it has been quite essential that the dynamic stochastic prob- 
lem has a Gibbsian limit, which physically is not totally unexpected due to 
the fluctuation-dissipation theorem. 

Having been interested mainly in proving the factorial form of the large-order 
asymptotes in the dynamic model considered we have not dwelled on the 
problem of effects of renormalization. However, in the static instanton analysis 
- according to [1,2] - the renormalization does non affect the solution of the 
stationarity equation and at leading order in 1/A^ it contributes to the C-type 
constants in correlation and response functions only. Details of the correct 
description of the renormalization in the instanton approach in the minimal- 
subtraction scheme can be found in Ref. [4]. Similar methods are needed for 
accurate calculation of the constants in our asymptotic expressions. 

As to the calculation of the dynamic critical exponent z, the renormalization- 
group equation for, e.g., the renormalized two-point function W2 = {^^)^ at 
the fixed point (7* 

W2 = 0, (35) 

may be used. In equation (35), k is the wave number and A is the scaling 
dimension known from the static theory. Details of the transformation of the 
g expansion considered above into the e expansion are described in Refs. [1,2]. 
From the factorial form of large orders in the e expansion the asymptotic 
estimate for the expansion of the exponent z follows: 

[^] kdM"^ + AMiyf ^ + Ai^^wl'^ - zmM""^ + ^^^^ ,,,, 

^ ~ m ) (36) 



a 


d 


1^- . - 


zt^ +A 


c)k 


dt 



where X^^'^ stands for corrections in 1/A^. Expression (36) demonstrates that 
the value z^^l has the same factorial form (34) as the static critical exponents. 
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